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Abstract 

We discuss various models and Chiral Perturbation Theory results for 
the Ki^ form factors F and G. We check in how much a simple parametriza- 
tion with a few parameters can be used to extract information from exper- 
iment. 
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1 Introduction 

In this note we give a simple parametrization and behaviour with few kinematical 
variables of the form factors contributing to the decay K + — > 7i + 7i^e + u. 

The interest of this work resides in the new measurements of K — > ttttCv (Kii) 
to be obtained with new facilities in the inmediate future. The increase in the 
number of events obtained will give more information about the dependence with 
the kinematical variables, generalizing and improving the previous and old results 
[[[]. This decay will be a major test of the consistency of the Chiral Perturbation 
Theory (CHPT), comparing the prediction of some 0(p 4 ) constants of the chiral 
expansion with the values given by other processes. Even more, this decay is 
expected to be the best one for a precise determination of the isoscalar tttt S- 
wave phase shift near threshold which will allow comparing CHPT and other 
approaches to chiral symmetry breaking. A review of the CHPT results can be 
found in || and references therein. 

One of the main results from the previous experiment is the observation of 
dependence on the energy of the dipion system (s n ). It was consistent with a 
linear dependence for all the form factors and with the same slope. The con- 
tribution from other kinematical variables was not considered. However, as we 
will see later, the expected variations of the form factors on the energy of the 
leptonic pair (si) and on the cos8 n , is on the same level as the errors in the 
previous experiment and will probably become important for the next generation 
of experiments. In |1] the form factors were fitted with a linear dependence only 
and the same slope for all form factors. The purpose of this note is to examine 
various theoretical expectations and check how well a simple extension of this 
parametrization performs. This will allow then an improved extraction of the tttt 
phase shifts from the data. 

Papers using CHPT and extensions for K^ are || and || . A review of early 
work is H . An alternative method to extract 7T7t phases is described in |J |7[] and 
references therein but it does not allow to use the full experimental information 
as much as our parametrization allows. 

The paper is divided as follows: we introduce the form factors in Sect. |[ 
We then discuss the singularities that appear in the various cuts in Sect. §. In 
Sect. ^ we study two different models for inclusion of resonances, in Sect. [5] 
CHPT at one loop and with inclusion of part of the p 6 effects in Sect. |[ We 
then discuss also some unitarization effects and finally the H form factor as well. 
Sect. H contains our proposal for the parametrization of the form factors. The 
last section contains our conclusions. 



2 Form Factors 

The decay K + (p) — > 7r + (pi)7r _ (p2)e + (^)z/(p !y ) is given by the amplitude 

T = ^§^"^^(1 - Tfe)^)^" - A") (1) 

where V M and A M are parametrized in terms of four form factors: 

TT 

A, = -i^-[P,F + Q,G + L,R} (2) 

with 

P = Pi + Pi , Q = Pi - Pi , L = pi + p u , iV = p e - p u . (3) 

and €0123 — 1-u and v are the lepton Dirac spinors and Gf is the Fermi constant. 

In this work we consider the theoretical predictions for the dependence of 
these form factors F, G and H, with the kinematical variables s n = (p\ + P2) 2 , 
s £ — ipe + Pu) 2 and cos^tt (9 n is the angle of the 7r + in the 7T7t system respect 
to the dipion line of flight). The .R-form factor is negligible in decays with an 
electron in the final state and we will not discuss it. 

The relation between form factors and the width can be found in several 
references 0, |3|, [(J . 

The main objective of this paper is to check in the existing models and cal- 
culations in how much simple parametrizations of F, G and H are sufficient. We 
will basically check if linear dependence on s n , s# and cos^ is sufficient or not. 

3 Cut Structure 

One indication that a constant or a linear structure is enough to describe a given 
form factor is, how far the masses of the possible intermediate states are from 
the physical domain. This is done by looking at the states contributing to the 
various cuts. The possible cuts with strong interaction intermediate states are 
shown in Fig. |l[ We discuss now the cuts in order of the figure. 

3.1 Cut (a) 

The strong interaction intermediate states that can couple here are in order of 
their mass: K, Kir, K*, Ktttt, K , K\ r . . . 

The K and Kq intermediate states can only contribute to the form factor R 
which is negligible in this decay. The Kn and K* intermediate state are vector 
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Figure 1: The possible cuts in the decay K + 
cut. The wiggly line is the weak current. 



7r + 7r e + v. The double line is the 



states and only contribute to the form factor H. Their effect has been estimated 
for Kir by one-loop CHPT and for the K* by a VMD estimate using the hidden 
gauge vector model in || . There was in effect a large cancellation between the two 
effects leading to a rather constant form factor across the relevant phasespace. A 
linear aproximation for both amplitudes separately should in any case be sufficient 
because the relevant momentum, sg is far from the thresholds. 

The K7rir has a threshold of about 0.765 GeV which is rather far from the 
physical region which ends at about 0.225 GeV. In addition in CHPT this only 
starts at the p 6 level. There are two K\ resonances, one at 1.27 GeV and one 
at 1.4 GeV. Both are again rather far from the physical region so we expect 
that their effect will also be small and certainly describable by at most a linear 
function of S£. 

If one checks the contribution of the axials of || we also notice that the 
axial- vectors only start contributing to this process at order p 6 . 

In conclusion, we expect the effects of this cut for K e4 decays to be linear and 
fairly small. 

3.2 Cut (b) 

Here the possible intermediate states are in order of mass again: im, p, /o, 

The 7T7T intermediate state is in fact one of the major reasons to study this decay 
in order to learn more about tttt phase shifts. The cut here is in the physical 
region for this decay. The effects of tttt final state rescattering will be discussed 
later, see also [§, [|. 

The effects of the p are confined mainly to the form factor G. There is a 
very small curvature in G due to this intermediate state visible in estimates of its 
effect, although in the region with sufficient data it is rather small. See section f| 
for its effect. 

3.3 Cut (c) 

This channel has charge 2 and has thus no resonance enhancement. The relevant 
intermediate state is mainly Kir and this also has a small phase shift in the 
relevant channel. We thus expect it to be well described by oneloop CHPT. 

3.4 Cut (d) 

Possible intermediate states here are: Kit, K*, Kq,. . . . The relevant variable in 
this channel is t and it is fairly far away from the threshold. Its effects are thus 
expected to be well described by a linear function in t, or linear in sg and cos8 n . 
Nonlinear effects on si are at most of order sj/m^^ which is about 1%. Actual 
estimates from CHPT at one-loop and various models confirm this expectation. 
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Figure 2: The form factor F and G of the model of Ref. [jTQ] as a function of s^ 
at cos 9 n = se = 0. 



3.5 Other singularities 

In addition to the above real singularities there are also kinematical singularities 
present in the various quantities. An example of this is the singularities present 
in the partial waves for Res,,- < 0. These follow because the values of t can 
become larger than tuk + fn n for some values of cos 8 n and sg. 



4 Resonances 



4.1 Resonances in a vector— axial- vector dominance model 



In this section we use the model of Ref. JTU| for r decays. Basically this model 
takes vector and axial-vector mesons and writes the amplitude such that in the 
low-energy limit all constraints from chiral symmetry are satisfied. In all channels 
contributions from the two lowest channels are included and the meson propaga- 
tors are described by Breit-Wigner functions including a s-dependent width. 

In Fig. H] we have shown the F and G form factors of this model as a function 
of s n for cos^ = se = 0. In accordance with the discussion in Sect. |3| we see 
that there is very little curvature. The curvature visible in G is due to the p-pole 
in this model. We have checked this by setting m p = 1.37 GeV in the model and 
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Figure 3: The form factor F and G of the model of Ref. [jTOj as a function of si 
at cos9n = and for J~s^ = 0.3 GeV and J~s^ = 0.35 GeV. 



then all curvature disappears. Experimentally there are very few events^ above 
s„ = 0.16 GeV 2 or E nn < 0.4 GeV. In the figure the fits linear in s n to the model 
curves below this energy are also shown. As can be seen the linear approximation 
is in the relevant region quite sufficient as an approximation to this model. 

Let us now turn to the Se dependence within this model. In Fig. |3] the 
form factors F and G are plotted for two values of s n . These values are in the 
region where accurate data can be expected. The lines are only plotted for the 
possible kinematical domains accessible in K e ±. As is obvious from the figure the 
S£ dependence in this model is very linear and any curvature can be neglected 
within the expected experimental accuracy. The differences in the slope for the 
two values of s n is small. 

We now turn to the cos^ dependence within this model. In Fig. § the form 
factors F and G are plotted for two values of s^. For G it is extremely linear 
while for F a very small curvature is present. The latter is however below the 
expected experimental uncertainties and show that the effect of -D-wavesQ is very 
small in this model. A linear best fit to F is also shown. The slope changes 



1 This will also be true in the next generation of experiments. We have in Fig. [lj shown the 
data of |jj] with the points plotted at the average energy in the bins used. 
2 Meant is here and below, D waves for F and higher than D waves for G. 
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Figure 4: The form factor F and G of the model of Ref. [1C] as a function of 

cos 6^ at se = and for J~s^ = 0.3 GeV and J~s^ = 0.35 GeV. 



somewhat between the two values of y/s^. The linear best fits are : 

F( v ^=0.3GeV) = 4.319 -0.191 cos 6 W 

F(^i; = 0.35 GeV) = 4.213 - 0.236 cos 6 W 

(2(^=0.3 GeV) = 4.429 -0.064 cos W 

G(y/s^ = 0.35 GeV) = 4.616 - 0.079 cos 6* (4) 

The differences in the cos 6^ part at the different energies are visible but do 
probably remain within the errors of the next generation of experiments. This 
can be seen in Fig. |] where we have for F also plotted the linear fit with the slope 
of the curve at ^ t /s^ = 0.35 GeV with the central value of the curve at ^/s^ = 
0.3 GeV. See also below for an approximate description of the s n dependence. 

So for this model a fit with linear s n , a linear cos 8 n and a linear si dependence 
seems sufficient in the region where there will be accurate data. 



4.2 Resonances using the model of Ref. 

Here we use the resonance parametrization as used in Ref. [^]. The difference 
with the previous section are that the resonances here only contribute at next-to- 
leading order in CHPT. We have therefore added the lowest order contribution 
m^:/(-i/(2)i 71 7r ) to F and G in order to make the plots more comparable to the 
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Figure 5: The form factor F and G of the model of Ref. 

at cos 6^ = Si = 0. 



as a function of s^ 



previous section. In this parametrization there is also no contribution from the 
axial-vectors to K ei , differences due to SU(3) breaking in the meson masses and 
couplings are also neglected but here scalars are included as well. In most cases 
however the main contribution comes from the vectors. 

The conclusions are also the same as in the previous section even if there are 
some differences in the numerical estimates of the various coefficients involved. 

In Fig. |5] we have shown the F and G form factors of this model as a function 
of s w for cos9 n = si = 0. In accordance with the discussion in Sect. |3| we see 
that there is again little curvature. The curvature visible in G is again due to 
the p-pole. We have checked this by setting m p = 1.37 GeV in the model and 
then all curvature disappears. Experimentally there are very few events above 
s n = 0.16 GeV 2 or E nn < 0.4 GeV. In the figure the fits linear in s n to the model 
curves below this energy are also shown. As can be seen the linear approximation 
is in the relevant region quite sufficient as an approximation to this model. 

In Fig. H we plot the sg dependence of the form factors F and G for two values 
of s n . The lines are only plotted for the possible kinematical domains accessible 
in K e4 . The S? dependence is very linear and any curvature can be neglected 
within the expected experimental accuracy. The differences in the slope for the 
two values of s n is small. 

We now turn to the cos8 n dependence within this model. In Fig. |7| the form 
factors F and G are plotted for two values of s n . For G it is extremely linear 
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Figure 6: The form factor F and G of the model of Ref. as a function of S£ 
at cos 9 n = and for ^fs^ = 0.3 GeV and ^Js^ = 0.35 GeV. 
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Figure 7: The form factor F and G of the model of Ref. as a function of 
cos 9 n at se = and for ^/s^ = 0.3 GeV and t/s^ = 0.35 GeV. 



while for F a very small curvature is present. The latter is however below the 
expected experimental uncertainties and show that the effect of .D-waves is very 
small in this model. A linear best fit to F is also shown. The slope changes 
somewhat between the two values of \fsZ. The linear best fits are : 



F(v^=0.3GeV) = 4.045 - 0.331 cos 6 n 

F(yfa = 0.35 GeV) = 3.940 - 0.399 cos 9 n 

G( v ^=0.3GeV) = 5.991 -0.230 cos n 

G(y/s^ = 0.35 GeV) = 6.195 - 0.281 cos 9 n (5) 

The differences in the cos^ part at the different energies are visible but do 
probably remain within the errors of the next generation of experiments. 

So for this model again a fit with linear s w , a linear cos^ and a linear se 
dependence seeme to be sufficient in the region where there will be accurate 
data. 



5 One-loop CHPT 

In Sect. |] we discussed in two models the effects of the poles present in the 
various cuts. Now we switch to a discussion of the effects of the continuum. 
In this section we will use CHPT to order p 4 . The contributions from the p 
constants Li to F and G are linear in the kinematical variables s n , t and u thus 
all curvature present is from the loop diagrams. To one-loop these diagrams 
have possible two-particle cuts, the main ones are the nir and the Kir ones as 
described in Sect. |[ Other cuts present in the loops are KK and Kr\. Their 
effect is however very linear in the kinematical variables. We use here the oneloop 
calculation of HI in the notation of H . For definiteness we use the values of the 



L\ as given in the DA$NE report ITTJ and all plots plot the real part only. Notice 
that this corresponds to the unitarized fit of Ref. so the fact that the plots do 
not agree well with the data of [[ij is due to the unitarization effects. 

In Fig. |8] we have shown the F and G form factors of this model as a function 
of s n for cos 6^ = se = 0. In accordance with the discussion in Sect. |3| we see 
that there is some curvature mainly generated by the tnr intermediate states. 
In the figure the fits linear in s n to the model curves below s n = 0.16 GeV 2 
are also shown. As can be seen the linear approximation in the relevant region is 
sufficient within the precision of the present experiment [[IJ . It becomes somewhat 
borderline for the next generation of experiments but the situation improves 
somewhat if we fit \F\ instead of the real part of F. The effect of the latter we 
have checked in the unitarized case where the curvature is more pronounced. 

Let us now turn to the si dependence. In Fig. |S| the form factors F and G are 
plotted for two values of s n . The lines are only plotted for the possible kinematical 
domains accessible in K ei . As is obvious from the figure the sg dependence in 
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Figure 8: The form factor F and G in p 4 CHPT as a function of s n at cos^ 
se = 0. 
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Figure 9: The form factor F and G in p 4 CHPT as a function of S£ at cos9 n = 
and for ^/s^ = 0.3 GeV and y^ = 0.35 GeV. 
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Figure 10: The form factor F and G in p 4 CHPT as a function of cosf^ at 
se = and for J~s^ = 0.3 GeV and J~s^ = 0.35 GeV. 



this model is again very linear and any curvature can be neglected within the 
expected experimental accuracy. The differences in the slope for the two values 
of Sir is small. 

We now turn to the cos9 n dependence within p 4 CHPT. In Fig. [10] the form 
factors F and G are plotted for two values of s n . For G it is extremely linear 
while for F a very small curvature is present. The latter is however below the 
expected experimental uncertainties and show that the effect of .D-waves is very 
small. A linear best fit to F is also shown. The slope changes somewhat between 
the two values of -Js^. The linear best fits are : 



5.138- 0.173 cos 0„ 
5.195 -0.225 cos ff 
4.858 -0.141 cos ff 
5.011- 0.177 cos^ 



F(v^ = 0.3 GeV) 
F(^ = 0.35 GeV) 

G(^=0.3 GeV) 
G{yfa = 0.35 GeV) 



(6) 



The differences in the cos (9^ part at the different energies are visible but do 
probably remain within the errors of the next generation of experiments. 

At all relevant values of s n the form factors are extremely linear in se and 
cos6 n . Parametrizing F as F = f s + fpCosO^ we get from one-loop CHPT 
\f p /f s \ = 0.043, 0.034, 0.008 at y/s^ = 350, 300, 280 MeV. This variation does 
probably stay within the future experimental errors. The curvature in s n might 
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Figure 11: The form factor F and G in partial p 6 CHPT as a function of s^ at 

cos 6^ = se = 0. 

become relevant in future experiments. 



6 Double Logs 



In this section we use the partial p 6 calculation of \T2\ for the F and G form 
factors. Here three types of contribution have been calculated using general 
renormalization methods. The relevant contributions are the ones proportional 
to (log(m|f / (J?)) 2 , Li \og(rriK //f) and LiLj. The scale \i is the subtraction scale 
in CHPT and we have chosen the scale of the logarithms to be the Kaon mass. 
The reason for this is that the same approximation in the p 4 calculation gives a 
reasonable agreement with the full result for this choice. Figures [II] to [13] show 
the same plots as in the previous sections. We have plotted here the partial p 6 
results with the lowest order added, the p 4 contribution is not included. The 
conclusions are basically the same as in the previous sections but notice that the 
curvature in F and G from this source is rather small. The contribution to the 
linear dependence on s n , S£ and cos 6^ is typically smaller than the p 4 results but 
not negligible. Absolute numerical results from this section are quite dependent 
on the particular choices of input [F2|. 
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Figure 12: The form factor F and G in partial p & CHPT as a function of Se at 
cos 6^ = and for ^/s^ = 0.3 GeV and y^ = 0.35 GeV. 
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Figure 13: The form factor F and G in partial p 6 CHPT as a function of cos^ 
at se = and for v /s^ = 0.3 GeV and ^J~s^ = 0.35 GeV. 
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7 Unitarization of F 

As we have seen in the plots, the behaviour of the form factors are almost linear in 
the models. Also, as we expected, the contribution from the S-w&ve from the im 
scattering produces some curvature in the F form factor that becames more clear 
in the case of the one- loop ChPT result. We will now discuss the S- wave part 
of this form factor, which is the source of the curvature using the unitarization 
method [[il| 14| . It consists in writing a dispersive relation for the partial wave 



using the ChPT limit to fix the subtraction constants. A previous study of the 
analytic properties separating the regions s n < and s n > provided a good way 
to introduce a general parametrization for the S- wave ||15|| , and the contribution 
from resonances in the t-channel (see || for details). 

The result is shown in Fig. [TJ]. The curve |/|jvii represents the absolute value 
for the S'-wave, where the contribution from the resonance sector, and therefore 
from the s n < region, is null. The other curves are the real and absolute value 
for the case with resonances. We have shown here both the absolute value and 
the real part for one case to show that the absolute value is more linear. The 
latter is the relevant one for fitting to the experimental data. 

A linear fit to the relevant energies for s„ is also plotted. The importance of the 
rescattering effects in the S"-wave is obvious, not only because of the modification 
of the linear behaviour, that can be explained by the fact that the cut is in the 
physical region, but also because of the contribution from the imaginary part. 

8 The H form factor 

The previous sections were devoted to the form factors F and G related to the 
cuts (b), (c) and (d) of the Fig. |l|. It remains to discuss the behaviour related 
to the H form factor. Using the 1-loop ChPT result, which is here order p 6 , 
and the hypothesis of resonance saturation for the unknown 0(p 6 ) constants, 
the behaviour of this form factor is not only linear but constant around the 
experimental value H{s lx = 0) = —2.68 ± 0.68 [1]]. There is a cancellation among 
the different contributions such that in the final curve the dependence with the 
energy is not appreciable, and the correction to the leading order is small M. 
Looking at the dependence with si and cos 9 n the shape is basically not modified 
with respect to the constant value. 

9 A sufficient parametrization 

In all the models/approximations analyzed in the previous sections we observed 
that F and G were always very linear in se with a slope that was very similar 
for both values of s^. We also observed that F and G were very linear in cos^ 
but with slopes somewhat varying between the two energies. The form factors 
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Figure 14: The S-wave in F using unitarization with (/#) and without (/nr) 
resonances as a function of s n at se = 0. The linearized fits for ^fs^ < 0.16 GeV 
are also plotted. 



depend on s n , t and u and the effects of the singularities other than the im ones 
can be expanded in t and u. By expanding to first order in these quantities we 
then indeed observe that both F and G have the structure 



W = w(s„) + w\Se + W2&-kX cos9 n . 
Where we used t — u = —2 <j v X cos 9 n with 



(7) 



cr„ 



1 1 - M , x = ^A 1 / 2 ^, S?r , s<) , A(o, 6, c) = (a - 6 - c) 2 - 46c . (8) 



Checking Eqs. (fD,© and @ we see that the s n dependence of the coefficient 
of cos 6* is indeed well described by <r n X. 

Using the partial wave expansions and their respective phases Sj (I: isospin, 
£: angular momentum), we then obtain as a parametrization f\ 



F = (fs(s«) + fis t )e<^ + f p a ir X cosO^^ 

G = (3p(s,)+9iSf)e^ w +3 rf (T T Icos^ 



(9) 



5 / p is a combination of the standard f p and g p partial waves 
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Then the curvature in s n in G is mainly due to the p resonance and in the 
experimentally relevant region is rather well described by a linear function in s n . 
We therefore advocate the use of 

G=(g P + g' s w + g t s e ) e^ 5 *' + g d a n X cos 9^^ (10) 



The phase S\ is relatively small and the last term in G is already a small cor- 
rection so one can neglect the difference between 5{ and 5® as well. For ^(s^) a 
linear approximation will be somewhat borderline, inclusion of a quadratic term 
is therefore useful and its presence should be checked in the data: 

/.(«*) = /. + /X + /X- (11) 

For H a linear approximation in s n should be sufficient 

H=(h p + h' pSn )e iSl ^K (12) 

So within this approximation we have 11 parameters plus the number of bins in 
s n for (<5q — Sl)(s w ) as fitting parameters. 

10 Conclusions 

We presented the results from several models for the form factors F and G 
and of CHPT calculations and extensions. We have shown that a simple lin- 
ear parametrization of the dependence on the kinematical variables is sufficient 
in most of the region where there will be sufficient data with the possible excep- 
tion of the ,S-wave part of the F form-factor. Even there a linear fit is quite a 
good approximation but then care has to be taken to define the theoretical slope. 
We showed that 11 parameters in addition to the phases are quite sufficient to 
fit the expected accuracy of the experiment and this should allow to determine 
those phases more accurately as a function of s w compared to the case where 
everything is fitted on a bin-by-bin basis. 
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